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I. INTRODUCTION 

Many processes in physics involve boundary surfaces which requires the solution of boundary 
and initial value problems. The introduction of a moving boundary into the physics usually pre- 
cludes the achievement of an exact analytic solution of the problem and recourse to approximation 
methods is required Q,Q] (see also P| and references therein). In the case of a moving plane bound- 
ary a time-dependent translation of the embedding space immobilizes the boundary at the expense 
of the increased complexity of the differential equation. It is the aim of this work to present an 
example of a soluble moving boundary and initial value problem in spherical geometry. 

The problems with moving boundary arise in many area of physics. One important example 
is sudden expansion of hot plasma with a sharp boundary in an external magnetic field which 
is particularly of interest for many astrophysical and laboratory applications (see, e.g., [J] and 
references therein). Such kind of processes arise during the dynamics of solar flares and flow of the 
solar wind around the earth's magnetosphere, in active experiments with plasma clouds in space, 
and in the course of interpreting a number of astrophysical observations [3-9] . Researches on this 
problem are of considerable interest in connection with experiments on controlled thermonuclear 
fusion (a recent review j^J summarizes research in this area over the past four decades). 

To study the radial dynamics and evolution of the initially spherical plasma cloud both analytical 
and numerical approaches were developed (see, e.g., [3-9] and references therein). The plasma 
cloud is shielded from the penetration of the external magnetic field by means of surface currents 
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circulating inside the thin layer on the plasma boundary. Ponderomotive forces resulting from 
interaction of these currents with the magnetic field would act on the plasma surface as if there 
were magnetic pressure applied from outside. After some period of accelerated motion, plasma 
gets decelerated as a result of this external magnetic pressure acting inward. The plasma has 
been considered as a highly conducting matter with zero magnetic field inside. From the point of 
view of electrodynamics it is similar to the expansion of a superconducting sphere in a magnetic 
field. An exact analytic solution for a uniformly expanding, superconducting plasma sphere in an 
external uniform and constant magnetic field has been obtained in The nonrelativistic limit 

n 

of this theory has been used by Raizer [lj] to analyze the energy balance (energy radiation and 
transformation) during the plasma expansion. The similar problem has been considered in L $] for a 
plasma layer. In the present paper we study the uniform expansion of the superconducting plasma 
sphere in the presence of a dipole magnetic field. For this geometry we found an exact analytical 
solution which can be used in analyzing the recent experimental and simulation data (see [10( and 
references therein). 

II. MAGNETOSTATIC TREATMENT 

In this section we first consider the simpler example of a non-relativistic expansion of the plasma 
sphere (v <C c, where v is the radial velocity of the sphere) in the presence of a dipole magnetic 
field. Consider the magnetic dipole p and a superconducting sphere with radius R located at the 
origin of the coordinate system. The dipole is placed in the position ro from the center of sphere 
(R < ro). The orientation of the dipole is given by the angle 9 P between the vectors p and ro- 
Here it is convenient to introduce the scalar magnetic potential ■i/'o ( r ) of the dipole magnetic field 
which is given by 

p. (r-r ) , * 

ipo(r) = - -3-. (1) 
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The dipole magnetic field is then calculated as Ho (r) = — V^o ( r )j 
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H (r) = -3 
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When the superconducting sphere is introduced into a background magnetic field the plasma 
expands and excludes the background magnetic field to form a magnetic cavity. The magnetic 



energy of the dipole in the excluded volume, i.e., in the volume of the superconducting sphere is 
calculated as 



Qr = ! 4^ dr ( 3 ) 
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where £ = R/vq < 1. This energy increases with decreasing 8 P and reach its maximum value at 
9 p = or P = it that is the magnetic moment p is parallel or antiparallel to the symmetry axis 
ro- In addition the magnetic energy Qr decays rapidly with the distance ro and for large ro ^> R 
is given by 



^ = ^(3cos 2 p + l). (4) 

In the case when the dipole approaches to the surface of the sphere ro ~ R the magnetic field of 
the dipole becomes very large and tends to the infinity as 



Qr = ^!_1±^£. (5 ) 
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We turn now to solve the boundary problem and calculate the induced magnetic field which 
arises near surface of the sphere due to the dipole magnetic field. Since the sphere is superconduct- 
ing the magnetic field vanishes inside the sphere. In addition the normal component of the field H r 
vanishes on the surface of the sphere. To solve the boundary problem we introduce the spherical 
coordinate system with the z-axis along the vector ro and the azimuthal angle <f> is counted from 
the plane (xz-plane) containing the vectors ro and p. Hence, using the expressions (|A2[1 - (|A4[1 the 
scalar potential Q at r < ro can alternatively be represented by the sum of Legendre polynomials 
(see the Appendix A for details) 
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The total magnetic field which is a sum of Ho (r) and the induced magnetic field is obtained from 
equation V • H = 0. Introducing the scalar potential, H (r) = — V?/> (r), the last equation becomes 



V 2 tp (r) = 0, i.e., ip (r) satisfies the Laplace equation. We must solve this equation with H = at 
r < R and boundary condition 
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We look for the solution of the Laplace equation which in a spherical coordinate system and at 
r > R can be written as 



i> ( r ) = V'o (r) + 
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where oli and Pi are the arbitrary constants and should be obtained from the boundary condition 
(JJJ). The second term in equation (jBJ) is the induced magnetic field. From equations one 
finds 
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Substituting equation Q into equation (JHJ) and using the summation formula obtained in Ap- 
pendix A from (JHJ) we find 



/ / \ p • Ro , Q ■ R* , , , s 
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where 



tpQD (r) = 
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Here r* = £ 2 r , Rq = r - r , R* 
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The term V'QD ( r ) in equation ()1U|) can be interpreted as a magnetic field of point-like quadrupole 
with the "quadrupole moment" D a p (r) and located in the j;z-plane inside the sphere at the distance 
r* (r* = £i? < R) from the centre. At large distances this term behaves as V'QD ( r ) — x zD xz /r 5 
with the quadrupole moment D xz = ^£ 5 psin# p {D aa = D xy = D yz = and a = x,y,z). The 



induced electric field is calculated from the Maxwell's equation V x E = — However, if plasma 
radial velocity is small, v/c <C 1, the amplitude of electric field is small as well (of the order of 
^Hq (r)) and may be completely ignored. Below we consider two particular cases for the magnetic 
dipole orientation in the space. 

(i) The case 6 P = 0; it. In this case the magnetic dipole is parallel or antiparallel to the vector 
ro- Obviously due to the symmetry reason the magnetic field does not depend on 4> and = 0. 
The magnetic field component Hg = — (1/r) (dip/d9) induces the surface current on the sphere. 
The ponderomotive forces resulting from interaction of this current with the magnetic field acts on 
the sphere surface with a magnetic pressure which can be calculated as an energy density of the 
magnetic field 
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This pressure vanishes at 9 = 0, it and has its maximum at 
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The value of # max tends to zero when the dipole comes close to the sphere and shifts toward the 
larger values, 9 max ~ tt/2, when the dipole goes to the infinity. Therefore the layer near 9 ~ # ma x 
of the expanding sphere will be mainly deformed by the external magnetic pressure. This behavior 
is clearly seen in the particle- in-cell simulation |l4j ]. 

The total force is calculated as a surface integral of the magnetic pressure 
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This force behaves as !F\\ ~ l~ s with s = 6 and s = 4 at large and small distances between the 
dipole and the surface of sphere, respectively. 

(ii) The case 9 P = tt/2. In this case there are two components of the surface currents which are 
proportional to Hg and at r = R. The magnetic pressure is then given by 
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where 
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At large distances, ( <C 1, the magnetic pressure is maximum at (j) ~ 5 and ^ (in equatorial 
plane), and = 0, 7T. At small distances, 1-(<1, only the region of sphere with 9 ~ 1 — £ ~ 
will be strongly deformed. 

The total ponderomotive magnetic force acting on the sphere is calculated as 



Again as for 6 p = 0, tt the force J-± behaves as J-± ~ l~ s with s = 6 and s = 4 at large and small 
distances, respectively. However, comparing equations (fT5|) and (fl9|) we conclude that the total 
magnetic force at 8 p = ir/2 is smaller than for parallel or antiparallel orientation of the dipole. For 
instance, from equations (|15|) and (|19|) we obtain J^y ~ 4:J T ± and J-\\ ~ 2J-± at £ <C 1 and £ ~ 1, 
respectively. 



In this section we consider the moving boundary problem of the plasma sphere expansion in 
the vacuum. In this sense unlike the magnetostatic problem considered above it is convenient here 
to introduce the vector potential of the induced and dipole magnetic fields. Consider a spherical 
region of space containing a neutral infinitely conducting plasma which has expanded at t = to 
its present state from a point source located at the point r = 0. The external space at the point ro 
contains a magnetic dipole p. The magnetic field of this dipole is given by Ho = V x Ao, where 
the vector potential Aq is 




(19) 



III. ELECTRODYNAMIC TREATMENT 



A 




(20) 



3 



r — 



As the spherical plasma cloud expands it both perturbs the external magnetic field and generates 
an electric field. Within the spherical plasma region there is neither an electric field nor a magnetic 
field. We shall obtain an analytic solution of the electromagnetic field configuration. 

We consider practically interesting case when the vectors p and ro are parallel (or antiparallel) . 
The general solution for the arbitrary orientation of p will be considered in a separate paper. 
Within this geometry the problem is symmetric with respect to the axis ro which is chosen as the 
axial axis of the spherical coordinate system. Then there is only one nonvanishing component of 
A , A 0r = A oe = 0, and 



|r - r | r Q l=1 V r o/ 

where P" (x) is the generalized Legendre polynomials with v = 1. Here D\ (x) = x l at x ^ 1 and 
Di (x) = x^ 1 ^ 1 at x > 1 as defined in Appendix A. 

Since the external region is devoid of free charge density, a suitable gauge allows the electric and 
magnetic fields to be derived from the vector potential A. Having in mind the symmetry of the 
original dipole magnetic field it is sufficient to choose the vector potential in the form A r = Ag = 0, 

oo 

A v (r, 9, t) = A 0ip (r, 6) + £ Ai (r, t) (cos 6) , (22) 

l=i 

and the components of the electromagnetic field are given by 



H - 19A ^ FT - 9A * F - 19A * (9* 



and = E r = Eg = 0. The equation for A\ (r, t) is obtained from the Maxwell's equations 



d 2 A + 2dA _ + _ l&A = Q (24) 
dr 2 r dr r 2 c 2 dt 2 

This equation is to be solved in the external region r > R(t) subject to the boundary and initial 
conditions. Here R (t) is the plasma sphere radius at the time t. The initial conditions are at t = 



(25) 



The first initial condition states that the initial value of A v is that of a dipole magnetic field. The 
second initial condition states that there is no initial electric field. Boundary conditions should be 
imposed at the spherical surface r = R(t) and at infinity. Because of the finite propagation velocity 
of the perturbed electromagnetic field the magnetic field at infinity will remain undisturbed for 
all finite times. Further, no incoming wave-type solutions are permitted. Thus, for all finite times 
Ai (r, t) — » at r — > oo. The boundary condition at the expanding spherical surface is H r = 
which can be replaced by A v (R (t) ,9,t) = or, alternatively, 



A l (R(t),t) = -^D l (^). (26) 

The problem of solving equation (|24|) subject to the initial and boundary conditions will be 
accomplished by the Laplace transform theory. The Laplace transform A\ (r, A) of the function 
Ai (r, t) is introduced by 
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Ai{r,X)= / Ai(r,t) e - Xt dt (27) 
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with Re A > 0. An inverse transformation is established by 



1 rcr+tco ^ 

A(r,t) = — Ai (r, A) e xt d\. (28) 
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The real parameter a should be larger than ReAj, a > ReAj, where Aj are the poles of A\ (r, A). 

The differential equation for A\ (r, A) is found from equations (j2*l|) and (|2*H|) and the initial 
conditions in Ef 



d 2 Ai(r,X) , 2dAi(r,X) 



dr 2 



+ 



dr 



1(1 + 1) Y 

r 2 c 2 



Al (r, A) = 0. 
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Its solution may be written as 



A t (r,A) = 4 

where hi (z) and (z) are the Hankel spherical functions and ai (A), c\ (A) are arbitrary functions 
of A determined from the boundary conditions. Since fi?' (z) gives rise to incoming waves, we should 
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set q (A) = 0. The solution to equation (|24|) at r > R (t) now may be written in the form 



oo 

r o 1=1 

where bi (A) = a/ (— iX). 

The moving boundary condition in equation (j'26f) requires the satisfaction of 
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Since the sphere moves with a radial velocity u less than the velocity of light c, we have R < ct 
or t — R (t) /c > 0. Thus, the contour in the integral of equation 1)32(1 should be closed by an infinite 
semicircle in the lower half plane and the integral evaluated by the method of residues. 

Explicit evaluation of this integral equation (|32|) . may be accomplished in the special case of a 
uniform expansion. Choosing the simple model of constant radial velocity R (i) = vt and assuming 
that R (t) < ro equation (j3*2*|) yields (see Appendix B for details) 



,, m _ (-1)' ("/ro)' jg 1 ,,,, 

where /3 = v/c < 1. Here P^ (z) are the generalized Legendre functions with z > 1, fj, = I, and 
u = -l-l. 

The solution of equations and <|31|) may be obtained by inserting equation Q33j) into (|3"Tj) 
and evaluating the integral (see Appendix B for details). The complete solution may finally be 
written in the form at v t < r < ct 



A, (r ,M) = W)-^g(^)'^P/(cos, 



(34) 



A v (r, 6, t) = Ao v (r, 9) at r ^ ct and (r, 0, t) = at r ^ -ui. Here £ = r/ci < 1, and 



(s) = 2'Z! (z 2 - 1) ^ Pf'- 1 (*) = ^ (r 2 - 1) ' dr. (35) 

The electromagnetic field components can be evaluated according to equation (|23|) . From equa- 
tions (|2SI) and JSU) it can be easily checked that the boundary condition on the moving surface, 



E (R) = -\ [v x H (R)) (or E v (R) = -f3Hg (R)), is satisfied automatically. It may also be noted 
that this special case of the uniform expansion falls within the conical flow techniques, as indicated 
in 3] for the case of uniform magnetic field. From symmetry considerations one seeks a solution of 
the form Ai (r,t) = r v § (r/ct). Substitution into the differential equation (|24|) yields an explicitly 
solvable ordinary differential equation whose solution, upon application of the boundary conditions 
($ (1) = 0, $ (J3) = -p/r 1 ^ 2 ), is given by equation 

It should be noted that all above results are valid only for R{t) < r$ or t < tq/v. At the 
time t = ro/v the dipole will enter into the plasma sphere and hence will be completely shielded 
by the latter. Therefore at t > tq/v the total electromagnetic field vanishes and the radiation is 
interrupted. 

IV. ENERGY BALANCE 

Previously significant attention has been paid P, ^| to the question of what fraction of energy 
is emitted and lost in the form of electromagnetic pulse propagating outward of the expanding 
plasma. In this section we consider the energy balance during the plasma sphere expansion in the 
presence of the magnetic dipole. When the plasma sphere of the zero initial radius is created at 
t = and starts expanding, external magnetic field Ho is perturbed by the electromagnetic pulse, 
H (r, t) = H (r, t) — Ho (r), E (r, t), propagating outward with the speed of light. The tail of this 
pulse coincides with the moving plasma boundary r = R (t) while the leading edge is at r = ct. 
Ahead of the leading edge, the magnetic field is not perturbed and equals Ho (r) while the electric 
field is zero. 

Our starting point is the energy balance equation (Poynting equation) 

V . S = -j.E-|-^±^, (36) 

J dt 8?r ' v ' 

where S = [E x H] is the Poynting vector and j = j'^e^, (with \e v \ = 1) is the azimuthal surface 
current density. The energy radiated to infinity is measured as a Poynting vector integrated over 
time and over the surface S c of the sphere with radius r c < ro (control sphere) and the volume f2 c 
enclosing the plasma sphere (r c > R or ^ t < r c /v). Integrating over time and over the volume 
f2 c equation can be represented as 



W S (t) = Wj (t) + AWem (t) , 



(37) 



where 
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W s (t) = lixrl \ dt' / S r sm9d6, Wj (t) = - dt' j • Bdr. 
Jo Jo Jo Jn c 



(38) 



Here S r = —£^E v Hq is the radial component of the Poynting vector. Wem {t) and AWem (t) = 
Wem (0)— Wem (t) are the total electromagnetic energy and its change (with minus sign) in a volume 
f2 c , respectively. Wj (t) is the energy transferred from plasma sphere to electromagnetic field and 
is the mechanical work with minus sign performed by the plasma on the external electromagnetic 
pressure. At t = the electromagnetic fields are given by H (r, t) = Ho (r) and E (r, t) = 0. Hence 
Wem (0) is the energy of the dipole magnetic field in a volume Q c and can be calculated from 
equation © by replacing R by r c and setting sin# p = 0, 



W EM (0) 
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where u = < 1. Then the change of the electromagnetic energy AWem (t) in a volume O c can 



be evaluated as 
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AWem (t) = - I - +H ~ H ° dr = Q («) - / ^^-dr. (40) 
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In equation ()40p f2^. is the volume of the control sphere excluding the volume of the plasma sphere 
(we take into account that H (r, t) = E (r, i) = in a plasma sphere). Hence the total energy flux, 
Ws (t) given by equation (|3*%|) is calculated as a sum of the energy loss by plasma due to the external 
electromagnetic pressure and the decrease of the electromagnetic energy in a control volume Q c . 
For non-relativistic (f3 <C 1) expansion of a one-dimensional plasma slab and for uniform external 
magnetic field (Ho = const) Ws — 2Wj ~ 2 A Wem, i-e., approximately the half of the outgoing 
energy is gained from the plasma, while the other half is gained from the magnetic energy . In the 
case of non-relativistic expansion of highly-conducting spherical plasma in the uniform magnetic 
field the outgoing energy Ws is distributed between Wj and AWem according to Wj = 1.5<5o and 
AWem = 0.5Q with W s = 2Q , where Q = H^R 3 /6 is the magnetic energy escaped from the 
plasma volume Therefore in this case the released electromagnetic energy is mainly gained 
from the plasma. 

Consider now each energy component Ws (t), Wj (t) and AWem (t) separately. Ws (t) is calcu- 
lated from equation (|38jl. In the first expression of equation (|38j) the t'-integral must be performed 



at 



^ t' ^ t [t < ^ ) since at ^ t' < ?f the electromagnetic pulse does not reach to the control 



surface yet and S r (r c ) = 0. From equations (|2*3*|) . (|3*4*|) and (|3*%|) we obtain 



W S (t) = Q(u) + 
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where r\ = r c /ct < 1. In non-relativistic limit, /3 — > 0, using the asymptotic expression (see, e.g., 
15]) pi (z) = z 2l+l / (21 + 1) at z -> oo, from equation (jUJ we obtain 
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with re = R 2 /ror c . In equation (|42[) Q (k) represents the magnetic energy of the dipole field in a 
sphere having the radius R* = R 2 jr c < R and enclosed in the plasma sphere. 

Next, we calculate the energy loss Wj (t) by the plasma which is determined by the surface 
current density, j. From the symmetry reason it is clear that this current has only azimuthal 
component and is localized within thin spherical skin layer, R — 5 < r < R + 6 with 5^0, 
near plasma boundary. Therefore in equation (|38|) the volume £l c can be replaced by the volume 
Qg ~ R 2 5 which includes the space between the spheres with r = R — 5 and r = R + 5. The surface 
current density is calculated from the Maxwell's equation, j = (1/4-zr) (cV x H — Within the 
skin layer we take into account that E = — - [v x H] and H r (R) = 0. Then 
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where 7 -2 = 1 — f5 2 and Sr are the relativistic factor and the surface of the expanding plasma, 
respectively. Note that the moving boundary modifies the surface current which is now proportional 
to 7~ 2 In equation (|43j) the term with ^-^j— has been transformed to the surface integral using 
the fact that the boundary of the volume Q$ moves with a constant velocity v and the electrical field 
has a jump across the plasma surface. Equation Q43JI shows that the energy loss by the plasma per 
unit time is equal to the work performed by the plasma on the external electromagnetic pressure. 
This external pressure is formed by the difference between magnetic and electric pressures, i.e., 




t t 

FIG. 1: The ratios Ys (t) (solid lines) and Tj (t) (dashed lines) for four values of (3 as a function of t (in 
units of ro/c) calculated from expressions (|41|l and l|44(l with r c = 0.5ro. 

the induced electric field tends to decrease the force acting on the expanding plasma surface. The 
total energy loss by the plasma sphere is calculated as 



Wj(t) 
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where £ = R/vq. In a non-relativistic case equation (|44j) yields: 
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The change of the electromagnetic energy in a control sphere is calculated from equation (|4()|) . 
At R < r c < ct (the electromagnetic pulse fills the whole control sphere) we obtain 



AW EM (t) = Q (u) 
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Comparing equations (|1T|) . and we conclude that AWem (t)+Wj (t) = Ws (t) as predicted 
by the energy balance equation (|H7j) . The non-relativistic limit of equation (|4H|) can be evaluated 
from equations (j4"2|) and 1)45 jl using the relation AWem (i) = W5 (i) — Wj(t). As an example 
in figure 1 we show the results of model calculations for the ratios T$ (t) = Ws (t) /Qo (t) and 
Tj (t) = Wj (t) /Qo (t) as a function of time (r c /c ^ t < r c /v). Here Qo (t) = Q (£) is the dipole 
magnetic energy escaped from the plasma sphere. For the relativistic factor (3 we have chosen a 
wide range of values. We recall that at ^ t ^ r c /c, i.e. the electromagnetic pulse does not yet 
reach to the surface of the control sphere, Ws (t) = 0. Unlike the case with uniform magnetic field 
discussed above (see also 0, Q]) there are no simple relations between the energy components 
Wg(t), Wj (t) and Qo(t). However, at the initial stage {t <C r c /v) of non-relativistic expansion 
the dipole field at large distances can be treated as uniform and the energies Ws (t) and Wj (t) are 
close to the values 2Qo (t) and 1.5Qo (t) (see figure 1), respectively. For any (3 the ratio V j (t) is 
almost constant and may be approximated as Tj (t) ~ Tj (0) or alternatively Wj (t) ~ 1.5CQo (t), 
where C = 7~ 6 (1 - /3) -4 (1 + 2/5) _/ is some kinematic factor. For (3 ~ 1 this factor is very 
large and behaves as C ~ (8/9) (1 — /3) _1 S> 1. As expected the total energy flux, Ws (t), increases 
monotonically with t. At the final stage (t = r c /v) of relativistic expansion (with (3 ~ 1) Ws — Wj. 
Hence in this case the radiated energy Ws is mainly gained from the plasma sphere. 



V. CONCLUSION 



An exact solution of the uniform radial expansion of a neutral, infinitely conducting plasma 
sphere in the presence of a dipole magnetic field has been obtained. The electromagnetic fields are 
derived by using the appropriate boundary and initial conditions, equations ()25j) and ()26|) . It is 
shown that the electromagnetic fields are perturbed only within the domain extending from the 
surface of the expanding plasma sphere r = R = v t to the surface of the expanding information 
sphere r = ct. External to the sphere r = ct the magnetic field is not perturbed and is given by 
the dipole magnetic field. In the course of this study we have also considered the energy balance 
during the plasma sphere expansion. The model calculations show that the radiated energy is 
mainly gained from the plasma sphere. For relativistic expansion the ratio Ws/Wj is close to 
unity and the radiated energy is practically gained only from plasma sphere. 

We expect our theoretical findings to be useful in experimental investigations as well as in 
numerical simulations of the plasma expansion into ambient nonuniform magnetic field. One of the 
improvements of our model will be to include the effect of the deceleration of the plasma sphere as 



well as the derivation of the dynamical equation for the surface deformation. A study of this and 
other aspects will be reported elsewhere. 
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APPENDIX A: SUMS WITH LEGENDRE POLYNOMIALS 



Using the known relation 



where D[ (x) = x l at |x| ^ 1 and Di (x) = at |x| > 1 one can derive some sums with Legendre 

polynomials Pi (cos 6) which are used in the main text of the paper. The first relation is obtained 
from equation (|A1|) by taking the partial derivative of the function Fq (x,8): 



Here the prime indicates the derivative with respect to the argument. 

The second relation follows from equation ()A1|) if we take the partial derivative over 6: 



a ■ a 00 
--F (x, 6) = - — - = J2 Dl (*) P{ (cos 0) , (A3) 

90 (i + x 2 - 2xcosey /2 fr{ 

where P^ (cos#) are the generalized Legendre polynomials PY (cos#) with u = 1. 
The third sum is calculated as 



J- [xF (x, 9)] = 1 a ^ = £ ^ (X) + Dl Pl (C ° S 6) ' (A4) 

ux (1 + x z — zxcosv) ' l=0 

Consider now the sum 



OO . r, OO . 

F (*> °) = E TTT xl+lp ' (cos 6) = ~dd E TTT xl+1Pl (cos 9) ' (A5) 



where x < 1. It is easy to see that 



a 2 



dxdO 

Using equation ()A6|) we finally obtain 



F (x, i 



d 



x sin ( 



dx 



(1 + x 2 — 2x cos 0) 



3/2 ' 



(A6) 



F(x,6») = sin6» 



x 2 F$ (x, 



x 2 sin ( 



F 3 (t, 0) tdt 



1 — X cos ( 



(1 + x 2 - 2xcos6) 3/2 sine y (I + x 2 - 2xcos6) 1/2 , 
In equation (jA7|) we have used the initial condition F (0, 9) = 0. 



(A7) 



APPENDIX B: EVALUATION OF THE VECTOR POTENTIAL 



For evaluation of the i nteg ral equation (|32j) we consider the explicit expression for the spherical 



Hankel functions h^p (z) [lj] 



' 'i\ k (l + k)\ 1 



k=0 



2 J k\(l-k)\z 



I 7 k+l 



(Bl) 



and assume that bi (A) = Bi/X l+1 , where F>i does not depend on A. This choice of b\ (A) assures 
that Bi is constant (see below). Inserting equation (|B1|) and bi (A) into equation (|H2|) we obtain 



k=0 



(l + k)\ {1-/3 



k\(l-k)\ V Pr 



k+l 



Qk+i+i (r) 



where r = i (1 — f3) > and 



(B2) 



3n (r) 



1 

2~TTi 



ia+oo e -i\r dX 1 QT, 



A n+1 n! <9<?" 



1 



A-<? 



9=0 



(B3) 



Here Img < ex. The integral within the square brackets according to the Kochi's theorem and at r > 
is equal to —e~ iqT . Therefore S y, (t ) = — (— ir) n jn\. Inserting this function into equation (|B2|) we 
arrive at equation fl-MJI (see, e.g., 15]). The complete solution is obtained by inserting equation 



into equation (|3Tj) and evaluating the contour integral as it was done above. 
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